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Abstract. We present radio-interferometric observations of HCN J = 1 — > line emission from the carbon star R Scl, obtained 
with the interim 3-mm receivers of the Australia Telescope Compact Array. The emission is resolved into a central source with 
a Gaussian FWHM of ~1", which we identify as the present mass loss envelope. Using a simple photodissociation model and 
constraints from single-dish HCN spectra, we argue that the present mass-loss rate is low, ~ 2 x 10~ 7 M Q yr -1 , supporting 
the idea that R Scl had to experience a brief episode of intense mass loss in order to produce the detached CO shell at ~10" 
radius inferred from single-dish observations. Detailed radiative transfer modelling yields an abundance of HCN relative to 
H2, /hcn> of ~ 10~ 5 in the present-day wind. There appears to be a discrepancy between model results obtained with higher 
transition single-dish data included and those from the J = 1 — » interferometer data alone, in that the interferometer data 
suggest a smaller envelope size and larger HCN abundance than the single-dish data. The lack of HCN in the detached shell, 
/hcn S2x 10~ 7 , is consistent with the rapid photodissociation of HCN into CN as it expands away from the star. 
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1. Introduction 

In the late stages of asymptotic giant branch (AGB) evolution, 
material is ejected via a slow (5-30 km s ) stellar wind into 
a circumstellar envelope of gas and dust. This mass loss (on 
the order of 10 7 -10" 5 M Q yr 1 ) is a key process for enriching 
and replenishing the interstellar medium, yet the mechanisms 
by which it occurs are not understood in detail. In particular, 
there is evidence that strong variations in the mass loss rate, 
perhaps related to flashes of helium shell burni ng, can lead to 
the fo rmation of circumstellar detached shells ( Olo fsson et alJ 
Il990l) . Such shells can be inferred from their excess emission 
at long infrared wavelengths (due to a lack of hot dust close to 
the star, Wille ms & de Jongl 19881) and confirmed by a ring-like 
morphology in maps of dust or CO line emission. 

Detached shells have been obser ved around about a 
half-d ozen AGB stars (se e review bv | Wallerstein & K napp 
1998b . b ut only U Cam jLindavist et alJ 1 19961 1 19991) and 
TT Cyg (Olofsson et al. 1998, 2000) have been mapped with 
millimetre-wave interferometers. Both cases reveal a beauti- 
fully symmetric shell showing uniform expansion from the star. 
High-resolution molecular line observations offer the possibil- 
ity not only to confirm detachment, but also to study variations 
in mass loss on short (<100 yr) timescales, related to the shell 
thickness, and to look for departures from spherical symmetry. 
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As one of the initial targets for the 3-mm receiver sys- 
tems being installed on the Australia Telescope Compact Array 
(ATCA), we have imaged the HCN / = 1 — » emission from 
th e carbon star R Scul ptoris (R Scl). Previous observations 
by Olofsson et al. ( 1996) with the Swedish-ESO Submillimetre 
Telescope (SEST) have indicated a detached shell which is only 
marginally resolved (radius «10") in the CO J - 3 —> 2 line, 
suggesting a relatively young age for the shell (< 10 3 yr, using 
the observed sh ell expansion velocity of 16 .5 km s _1 and a dis- 
tance of 360 pc: ISchoier & Olofssonfcoo'll) . Detections of other 
molecules like HCN, CS, and CN have also been reported by 
lOlofsson et alJ Jl996i) . These molecules are more easily pho- 
todissociated than CO, but it is not clear whether they are lo- 
cated in the detached shell or confined to the current mass loss 
envelope, which the SEST data, with limited spatial resolution, 
could not distinguish from the shell. Another outstanding ques- 
tion is the relationship between the molecular gas shell and the 
larger dust shell (radius ~20") seen in scattered optical light 
dGonzalez Deleado et al J200lU2003l) . The recent commission- 
ing of a southern millimetre array motivated us to obtain much 
higher resolution data toward this object. 

In this paper we present the results of our ATCA observa- 
tions and analyse them, along with new and previously pub- 
lished single-dish spectra, using a detailed r adiative transfer 
code for modelling circumstellar line emission JLindqvist et alJ 
2000: ISchoier& Olofss onl200ll) . 
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Table 1. Single-dish HCN molecular line observations of R Scl. 



-100 



Fig. 1. Coverage of the visibility plane obtained after com- 
bining observations in all four ATCA configurations (EW214, 
EW352, H168, and 750A). 



2. Observations and data reduction 

2.1. Interferometer data 

At the time of the observations, the ATCA had three antennas of 
22 m diameter equipped with dual polarisation 3-mm receivers 
covering the bands 84.9-87.3 and 88.5-91 .3 GHz. We observed 
R Scl between 2002 June and October in four different array 
configurations: EW214 (baselines of 31, 61, and 92 m), EW352 
(46, 77, and 122 m), H168 (61, 111, and 171 m), and 750A 
(138, 275, and 413 m). Note that in the H168 configuration, the 
3-mm antennas are situated in a north-south line; the remain- 
ing configurations are east-west arrays. The 750A observations 
were taken over two nights spaced 3 days apart; observations 
in the other configurations occurred on single nights. The com- 
bined visibility plane coverage is shown in Fig. [2 

All observations were conducted in mostly clear weather, 
with above-atmosphere single-sideband system temperatures 
of T* » 300 K near the zenith. The correlator was config- 
ured to receive both linear polarisations in two frequency win- 
dows: a narrowband (spectral line) window centred on the HCN 
/=!_>() line at 88.6 GHz with 128 channels across 32 MHz, 
and a wideband (continuum) window centred at 89 GHz with 
32 channels across 128 MHz. The pointing and phase centre 
was at ff 2 ooo = 01 h 26 m 58 s .49, 5 2 ooo = -32°32'35'/5, which is 5" 
east of the stellar position (in the subsequent data analysis we 
have applied an offset so that the phase centre coincides with 
the star). At 89 GHz the ATCA primary beam has a FWHM of 
about 35". 

Gain calibration was performed by frequent observations 
(every 10-20 minutes) of the nearby quasars B0 104-408 or 
B0208-512. The amplitude and phase gains were derived us- 
ing the continuum window data and then transferred to the 
spectral line window after applying offsets determined from a 
half-hour integration on the bandpass calibrator (B1921-293, 
B2255-282, or B0420-014). Channel-dependent gains were 
also determined from the bandpass calibrator. Uranus was used 
to set the flux scale, assuming a uniform disk with brightness 
temperature of 134 K. We also adjusted the antenna pointing 
once an hour on B0208-512; typical pointing shifts were ~5". 
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a Total integrated intensity calculated over full extent of line. 
b Estimated from a Gaussian fit to the observed spectrum. 
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All data processing was conducted using the MIRIAD 
package. As an initial step in the reducti on, values of T* mea- 
sured via the chopper wheel method JKutner & Ulichlfl98lh 
were used to rescale the data, then an elevation-dependent gain 
curve was applied based on observations of an SiO maser taken 
on 2002 June 4. Subsequent processing was as described above 
using standard routines. The calibrated visibilities were Fourier 
transformed using robust weighting and a (X'5 pixel size with a 
channel spacing of 1 km s _1 (roughly the effective velocity res- 
olution given the original channel spacing of 0.25 MHz a 0.8 
km s _1 ). The maps were then cleaned down to a 2<x level over 
the inner 20" x 20". 

Over the course of the observations, it became clear that 
there were residual phase errors in the data that were not be- 
ing removed by the standard technique of phase referencing off 
a nearby quasar. These phase errors, revealed as phase jumps 
when switching sources, can be modelled as errors in the an- 
tenna positions of up to a few mm. For most of the data, we 
were able to minimise the effects of these phase errors by refin- 
ing the baseline solution from observations of multiple quasars 
taken on the same or adjacent days. No such correction could 
be applied to the EW352 data, which are therefore excluded 
from the image analysis below (Sect. 13 . 11 . but are included in 
the visibility analysis after self-calibration (Sect. 1331 . 

2.2. Single-dish data 

As a complement to the interferometer observations, multi- 
transition single-dish data have been obtained. The single- 
dish data are used to check that the interferometer recovers 
all of the source flux and also to constrain the circumstellar 
HCN models. Using the SEST, the HCN J = 1 -> and 
J — 3 — > 2 transitions were observed in February 2003 and 
December 1996 respectively. Two acousto-optical spectrome- 
ters were used at SEST (86 MHz bandwidth with 43 kHz chan- 
nel separation, and 1 GHz bandwidth with 0.7 MHz channel 
separation). Dual beam switching (beam throws of about 11'), 
in which the source was placed alternately in the two beams, 
was used to eliminate baseline ripples. The line intensities are 
given in the main beam brightness temperature scale (T m \,), i.e., 
the antenna temperature has been corrected for atmospheric at- 
tenuation (using the chopper wheel method) and divided by the 
main beam efficiency. In addition, re cently published HCN / = 
3 — > 2 and J - 4 —> 3 spectra from Bieging (2001), obtained 
with the Heinrich Hertz Submillimeter Telescope (HHT), have 



T. Wong et al.: ATCA imaging of HCN emission from R Scl 



3 



-35.5 



-25.5 



-15.5 



-5.5 



-33.5 



-23.5 



-13.5 °0 



O 




-3.5 



-31.5 







-21.5 



-11.5 



-1.5 



-29.5 



-19.5 



-9.5 



0.5 



■0 



-27.5 



t-17.5 o<2 



/ 



-7.5 



2.5 



10 -10 10 -10 10 -10 10 -10 10 -10 

RA offset (arcsec; J2000) RA offset (arcsec; J2000) RA offset (arcsec; J3000) RA offset (arcsec; J2000) RA offset (arcsec; J2000) 

Fig. 2. Deconvolved channel maps of HCN emission from R Scl. The contour levels are 0.05n 2 Jy beam 1 , where n = 1, 2, 3,4, 
and the beam size is 2'.'7 x l'/8. The velocity channels (given in the LSR frame) have been binned to 2 km s _1 . 



been used. The single-dish data are presented in Table ^ an d 
Sect. IP 

3. Observational results 

3.1. Channel maps 

Figure |2] shows the HCN J = 1 — > channel maps (with 
the EW352 data excluded) binned into 2 km s channels. 
The synthesised beam is 2'.'7 x 1'.'8 with a position angle of 
18°. The RMS noise in the 1 km s" 1 channel maps is <x = 23 
mJy beam -1 ; the peak emission reaches a level of about 40<r. 
We note here that the HCN J - 1 — > line actually consists of 
three hyperfine components of relative strength 3:5:1 (in order 
of increasing frequency), with the low-frequency satellite line 
shifted +4.8 km s and the high-frequency satellite line shifted 
-7.1 km s _1 with respect to the main line. 

Inspection of the channel maps suggests that the source is 
moderately resolved, a result which is confirmed by the visibil- 
ity analysis in Sect. 13.31 The source also appears to be some- 
what extended to the east. We believe that this extension is most 
likely attributable to the instrumental phase errors discussed in 
Sect. 12. II mapping the EW214 data alone shows the emission 
to be displaced by about / .'6-077 to the east of the star's posi- 
tion, whereas the other tracks do not show such an offset. The 



phase error needed to produce such an offset is ~30°, which is 
within the range expected from residual errors in the baseline 
solution. 

Fitting two-dimensional Gaussians to the individual chan- 
nel maps, we find no systematic variations in source size or po- 
sition angle with velocity. Across the peak of the line profile, 
from -10 to -20 km s~', the average values of the major and 
minor axis FWHM are 3'.'1 and 2'.'3 respectively, with a posi- 
tion angle of 36° (note that the beam has not been deconvolved 
from these values). 

We see no indication of an expanding shell structure corre- 
sponding to that seen in the SEST CO / = 3 — » 2 data. A shell 
of radius 10" would appear well within the half-power radius 
of the primary beam, especially on the eastern side, since the 
pointing centre was offset 5" east of the star. Moreover, an ex- 
panding spherical shell would become smaller near the extreme 
velocities (aside from some complications due to the hyperfine 
structure of the HCN line), and thus certainly fall within the 
primary beam, yet there is no indication that the source struc- 
ture or size changes with velocity. 

3.2. Spectra 

The integrated HCN ./ = 1 — > ATCA spectrum is shown 
in Fig. |3 a), overlaid on the equivalent SEST spectrum. A 
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Fig. 3. (a) HCN (1-0) spectra for R Scl, with the dashed 
line derived from the ATCA image (integrating the channel 
maps over a 6" x 6" box centred on a 2 ooo=01 h 26 m 58 s .09, 
5 2000 =-32°32'35'.'5) and the thin solid line from the SEST 
observations, (b) Average visibility phase spectrum, including 
data from all baselines and array configurations. The phase vs. 
time has been self-calibrated to bring the average phase to 0°. 



Gaussian fit to the ATCA spectrum gives a mean velocity, 
V*, of -15.8 km s~ l and a FWHM of 19.1 km s -1 ; how- 
ever, this may be affected by the hyperfine structure of the 
HCN line. The peak flux of 1.2 Jy agrees, within the calibra- 
tion uncertainties, with the single-dish flux of 1 .4 Jy detected 
with the SEST (S = 7] B T mb r-\ where ^=0.75, r mb =75 mK, 
and F~'=25 Jy Kr 1 ). Thus, our lack of interferometer spacings 
less than 30 m does not appear to have resulted in a significant 
amount of "missed" flux, at least within the 58" SEST beam. 
This confirms that virtually all of the HCN emission comes 
from the central component. 

The higher-transition spectra taken with the SEST and 
HHT (see Tabled are in general agreement with the / = 1 — » 
spectra, except that their mean velocity is » -17.5 km s _1 and 
the linewidths are somewhat narrower, FWHM « 15 km s . 
Since these lines are less affected by hyperfine structure, they 
should serve as a better indication of the kinematics of the en- 
velope. The spectra are presented and discussed in Sect. 14.21 
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Fig. 4. FWHM of axisymmetric Gaussian fits to the visibility 
data averaged in 3 km s _1 channels. The source is constrained 
to lie at the phase centre. The flux of the model Gaussian, on 
the same numerical scale but in Jy, is overplotted as a dashed 
line. 
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3.3. Visibility analysis 



The results of Sect.|3J]indicate that the HCN J = 1 -> emis- 
sion comes from a compact central source that is only moder- 
ately resolved by our observations. We now proceed to analyse 
the visibilities directly, under the assumption that the emission 
structure is axisymmetric, i.e. there is no phase information 
in the visibilities. This assumption allows us to self-calibrate 



UV Distance (kX) 

Fig. 5. (Top) Averaged visibility amplitudes, in 5 kA bins, 
across the peak of the line profile (filled circles). The solid line 
is the model fit, with a FWHM of l'.'O. The open squares rep- 
resent the averaged visibility amplitudes of the gain calibrator, 
B0104-408. (Bottom) Same but for the high-velocity edge of 
the line profile. The solid line is the model fit, with a FWHM of 
1'.'4, whereas the dotted line corresponds to a FWHM of l'.'O. 
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the data in phase (in practice a 30-minute interval was used) 
and therefore include the EW352 data, which were strongly af- 
fected by phase errors. We stress that applying self-calibration 
with only three interferometer baselines is guaranteed to pro- 
duce an axially symmetric source at the phase centre, as long as 
closure errors are negligible. Independent evidence that there is 
little additional source structure comes from the averaged visi- 
bility phase spectrum shown in Fig.|3jb). Since the channel-to- 
channel phase has not been self-calibrated, but only calibrated 
using the bandpass calibrator, the flatness of the phase spectrum 
indicates that the source structure is the same at all velocities. 

We fit a model of an axisymmetric Gaussian to the self- 
calibrated visibilities. The resulting FWHM in arcseconds is 
plotted in Fig.0] along with the flux of the model. Most of the 
values are close to l'.'O, although somewhat higher values (l'/3) 
are seen at the redshifted edge of the line. An average over all 
points gives 1'.'05 + 0717. Note that a FWHM of 1" corre- 
sponds to a radius at half maximum, at the distance of R Scl, 
of about 200 AU (3 x 10 15 cm). As with the image analysis in 
Sect. 13.11 we do not see the variation in size with line-of-sight 
velocity expected for an expanding shell. 

The top panel of Fig. |5] shows a 170 Gaussian model fit 
to the peak of the line profile, overlaid on the visibility ampli- 
tude averaged in 5 k/i bins. For comparison, the open squares 
give the average visibility amplitude for the phase calibrator, 
B0104-408, during the two October runs in 750A (the ampli- 
tude was ~3Q% lower in the June runs, but this may be due to 
intrinsic variability of this quasar and/or errors in relative flux 
calibration). Note that the amplitude of the calibrator does not 
show a significant decrease with baseline length, as would be 
expected if atmospheric phase decorrelation were affecting the 
averaged visibility amplitudes. Consequently, there is no doubt 
that the HCN emission from R Scl has been spatially resolved 
by our observations. However, because our baselines are not 
long enough to show a possible null in the visibility function, 
other types of extended models besides a Gaussian (e.g., a uni- 
form disk or optically thin sphere) are not excluded. 

The bottom panel of Fig.[5]shows a 1736 Gaussian model 
fit to the redshifted edge of the line profile, which showed in- 
dications of a larger source size in Fig.|4] While the fit appears 
to be a reasonable one, there are clearly some bins that agree 
better with the 1" source size derived from the full velocity 
range (dotted line). Thus, the significance of the larger source 
size at these velocities appears doubtful. We note that for an 
expanding shell that has been well-resolved, we would expect 
the source size to decrease, rather than increase, near the line 
wings. 

4. Radiative transfer modelling 

To model the circumstellar line emission we have used a de- 
tailed non-LTE radiative transfer code, based on the Monte 
Carlo m ethod fe ernes|| 1 9791) . The code has been describ ed in 
detail in Schoier (2000^1^ ISchoier & Olofssonl feOOll) and 
has been benchmarked a gainst a wide variety of mol ecular line 
radiative transfer codes in lvan Zadelhoffetal] ( l2002l) . The code 
assumes a spherical envelope expanding at constant velocity 
and includes radiative excitation through vibrationally excited 



states, and a full treatment of line overlaps between various hy- 
perfine components. Ou r procedure and the re levant molecular 
data are summarised in lLindqvist et alJ |2000). 

From the HCN line emission alone it is not possible to con- 
strain both the abundance of HCN (relative to H2) and the mass 
loss rate. For our analysis, therefore, we first adopt a mass loss 
rate for the present day wind. We construct a model of the CO 
envelope to derive the temperature profile, then determine the 
best-fit values for the HCN abundance and envelope size, us- 
ing both the single-dish and interferometer data as constraints. 
The HCN abundance, /hcn, is taken to be fixed throughout the 
envelope out to a specified outer radius, thcn- 

4.1. The mass loss rate and temperature profile 

lLeBertrel(ll997l) estimated the present mass loss rate for R Scl 
to be approximately 10~ 7 M yr , scaled to our adopted dis- 
tance of 360 pc, based on radiative transfer modelling of the 
observed dust emission. This method assumes that the dust 
grain properties and the dust-to-gas mass and velocity ratios 
in R Scl are the same as in IRC+10216, for which the mass 
loss rate is taken to be 2.5 x 10~ 5 M yr -1 at a distance of 
200 pc. A larger mass loss rate of 4 x 10~ 7 M n yr" 1 (a gain 
scaled to our adopted distance) is derived by Gustaf sson et alJ 
d!997l) . using Ki scattered light observations, but this esti- 
mate must be considered relatively uncertain since it is based 
on a method which has not been extensively tested. As long 
as the size of the emitting region is limited by photodisso- 
ciation rather than excitation, a rough estimate of the mass 
loss rate can also be obtained from the observed size of the 
HCN brightness distributio n, using a simple photodissocia- 
tion model (cf.. iLindavist et afll2000l) . For an estimated ex- 
pansion velocity of 10.5 km s _1 (see below), a mass loss rate 
of approximately 5 x 10~ 7 M yr~' is derived assumin g the 
same dust properties as for IRC+10216 (Sch oier & Olofssonl 
2001). Based on these estimates, we adopt a mass loss rate 
of 2 x 10~ 7 M yr~' for the present epoch. Thus, the mass 
loss rate for R Scl appears higher than the present mass loss 
rates of sources with rela tively old detached shel ls lik e TT Cyg 
(M ~ 3 x 10~ 8 M^ vr~' : loiofsson et al.lll998l EoOO) but com- 
parable to that of the 'young shell' source U Cam (M ~ 
2.5 x 10~ 7 M„ vr~': ILindavist et alll996lll999t) . 

The expansion velocity of the wind is a free parameter in 
the radiative transfer that is well constrained by the width of 
the observed lines. A value of 10.5 km s _1 can reproduce both 
the J = 1 — » line and the higher J transition HCN line 
profiles, which are much less affected by hyperfine structure 
(see Sect. 14.21 . The source is assumed to have a distance of 
360 pc and the stellar radiation to be described by a blackbody 
with a luminosity of 5500 L and a temperature of 2500 K 
(Schoier & Olofsson 2001, Kerschbaum priv. comm.). The in- 
ner radius of the envelope is fixed at 3 x 10 14 cm (20 AU). 

Using the above parameters we have computed a CO model 
of the envelope, adopting a CO abundance of 1 x 10~ 3 relative 
to H2. The outer radius of the CO envelope, ~3 x 10 16 cm (6"), 
is ob tained from a CO photodissociation model JMamon et alJ 
1988). In addition to determining the steady-state level popu- 
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Fig. 6. x 1 ma P showing the sensitivity of the HCN envelope 
model to the adjustable parameters ,/hcn and thcn- Contours 
are drawn at x^ in + (2.3,6.2, 11.8) indicating the 68% ('lcr'), 
95% ('2<t'), and 99.7% ('3<x') confidence levels, respectively. 
The quality of the best fit model (Model 1) can be estimated 
from the reduced chi-squared statistic Xr e d=XmiJ(N-2), and is 
shown in the upper right corner. Also shown is the number of 
observational constraints used, N, which includes the single- 
dish spectra and the ATCA spectrum at the centre pixel. The 
model that best reproduces the visibility amplitudes obtained 
by ATCA (Model 2) is also shown for comparison. 

lations of the CO molecule, the code simultaneously solves the 
energy balance equation, including the most relevant heating 
and cooling processes. Heating is dominated by collisions be- 
tween the dust and gas, while cooling is generally dominated by 
CO line emission, although adiabatic cooling due to the expan- 
sion of the envelope is also important. We have assumed that 
R Scl has the same dust properties as IRC+10216. Our excita- 
tion analysis allows for a self-consistent treatment of molecu- 
lar line cooling. In this way the kinetic temperature structure of 
the envelope is determined, for use in the subsequent HCN line 
modelling. Throughout most of the envelope the temperature 
structure is well described, to within 10%, by a simple power- 
law Tan = T (r /r) 0A7 , where T = 182Kandr = lxl0 16 cm. 
However, in the inner parts of the envelope (r < 1 X 10 15 cm) 
the temperature rises more dramatically to a value of 1350 K at 
the inner edge of the envelope. 

We find that our derived CO line intensities are much lower 
than the observed ones, consistent with most of the observed 
CO line emission coming from a detached shell distinct from 
the present mass loss envelope. A more thorough discussion 
will be presented in Lindqvist et al. (in prep.). 

4.2. HCN line modelling 

With the physical structure of the envelope fixed, a set of mod- 
els was calculated where the HCN abundance relative to H2, 



/hcn, was allowed to vary from 5 x 10~ 6 to 5 x 10~ 5 and the 
outer radius of the HCN-emitting envelope from 2 x 10 15 cm 
to 5 x 10 16 cm. For each m odel a total x 1 (f° r a definition 
see Schoi er & Olofssonl200ll) was determined via comparison 
with the observed spectra, assuming a calibration uncertainty of 
20% for each spectrum. In addition to the single-dish data pre- 
sented in Table [2 the ATCA spectrum at the centre pixel of the 
maps was used in the^- 2 analysis. Note that the^ 2 map (Fig.|6) 
shows that the line intensities become relatively insensitive to 
the adopted size of the HCN envelope when > 3 x 10 16 cm. This 
is where the emitting region becomes excitation limited. Such 
solutions, however, are implausible given that the CO envelope 
is calculated to extend only to 3 x 10 16 cm, and CO is less 
readily photodissociated than HCN. The best fit model has an 
abundance of (1.2+0.4) x 10~ 5 and a size of ( 1 .4+0.5) x 10 16 cm 
(Fig.|6j Model 1). The reduced^ 2 has a value of 1.3, indicating 
a good fit. However, as can be seen by comparing this model 
with the observed spectra (Fig. [7}, the model ATCA spectrum 
is too weak and presents more of a double-peaked rather than 
parabolic line shape. Note that we have not included the vibra- 
tionally excited lines seen in the J — 3 — > 2 and ./ = 4 — > 3 
spectra in Fig.0in our analysis, given that they are masing and 
highly variable (Bieging 2001). 

The main source of the discrepancy between our model and 
the ATCA data is that the size of the model envelope is approxi- 
mately five times larger than that obtained from the interferom- 
eter Mv-plane data (Sect. l3~3l . Applying the ATCA uv sampling 
to the model envelope, we find that the model severely under- 
predicts the flux observed on long baselines (Fig.|SJ. As shown 
by the dashed line, a better fit is obtained by raising the abun- 
dance of HCN (/hcn — 4.5 x 10~ 5 ) while making the envelope 
smaller (recN = 5.5 x 10 cm) (Model 2). However, such a 
model would predict much stronger emission (by ~50%) in the 
higher J transitions, in conflict with the single-dish constraints 
(Fig.|6j. We also note that Model 2 predicts a decrease in source 
size at the extreme velocities of the line profile, which, while 
expected for a spherically symmetric wind expanding at a con- 
stant velocity, is not evident in the observations (Fig.0}. 

It is not clear how to resolve the discrepancy between the 
modelling and data. Changing the adopted mass loss rate af- 
fects the HCN abundance but not, to a first approximation, the 
HCN envelope size. Since the HCN lines are mostly radiatively 
excited, mainly through the 3-/im stretching mode which lies 
close to the peak of the spectral energy distribution, chang- 
ing the kinetic temperature structure within reasonable limits 
also does not significantly affect the model results. While the 
HCN / = 1 — > hyperfine lines can be inverted in parts of the 
envelope, the effect of maser emission on the scales observed 
here (> 1") should be small. Previous mode lling of H 12 CN 
and H 13 CN data for a sample of carbon stars bv lLindavist et alj 
(2000) indicates that our technique may be less successful at 
modelling lines of high optical depth such as HCN, where 
changes in the physical parameters in the envelope, such as 
density structure and velocity field, will have large effects. 

A higher degree of sophistication in the modelling would 
require more detailed observations. The near-term prospect of 
higher resolution HCN / = 1 — > observations is remote, but 
higher resolution mapping of other molecules (such as CO and 
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Fig. 7. Best fit model (full line), using a HCN abundance of 1.2X 10~ 5 and envelope size of 1.4x 10 16 cm, overlaid on the observed 
spectra (histograms). The narrow spectral features in the J — 3 — > 2 and J = 4 — > 3 spectra are from vibrationally excited HCN. 
The ATCA spectrum is taken at the centre pixel of the image cube. 



CN) would be useful. In particular, good spatial resolution will 
be important in order to address the possibility that the medium 
is hig hly clumped, as has been suggested by Olo fsson et alJ 



5. Discussion: The detached shell 

Since our HCN data do not show direct evidence for a de- 
tached shell around R Scl, it is worth reviewing the evidence 
that such a shell exists. First there are the far-infrared colours 
of this star, which are similar to those of stars with clearly de- 
tached envelopes (Olofsson et al. 1990). In addition, the CO 
J — 3 — > 2 SEST maps presented by Olofsso n et alJ (11996), 
with an effective resolution of about 7", show the systematic 
change in source size with velocity expected from an expand- 
ing spherical shell, as well as off-centre peaks in the central ve- 
locity cha nnels that are indicative o f detachment. Third, as dis- 

I | I I 

cussed by Gonzalez Delgado et al. (2003), attempts to model 

the CO data with only a continuous mass loss envelope lead 
to very high mass loss rates (~10~ 5 M Q yr _1 ), which cannot 
be reconciled with the double-peaked line profile nor produce 
the extended brightness distribution of the CO J = 3 —> 2 
line. A high mass loss rate is also inconsistent with the fact that 
the star is optically bright. Finally, the detection of a geometri- 
cally thin dust shell in po larised optical light is also noteworthy 
jGonzalez Delgado et alJl2003l) . although this shell appears to 
have a radius roughly twice as large as that of the proposed CO 
shell. 
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Fig. 8. Visibility amplitudes of the observed HCN J = 1 — > 
line emission obtained at ATCA towards R Scl as a function 
of the projected baseline length. The observations, integrated 
from -28 to -4 km s~' and binned to 5 Li, are plotted as filled 
symbols with Icr error bars. Also shown is the result of apply- 
ing the same uv sampling to the best-fit circumstellar model 
including the single-dish data (Model 1; solid line). A model 
with higher abundance and smaller envelope size is able to re- 
produce the interferometer observations much better (Model 2; 
dashed line). 



The modelling presented in this paper provides indirect 
support for the hypothesis that R Scl has experienced a brief 
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episode of intense mass loss, leading to the formation of a 
detached CO-emitting shell. Our derived HCN abundance of 
about 1 x 10~ 5 is in excellent agreement with the estimated 
photosph eric HCN abundance o f ~ 8 x 1(T 6 for R Scl deter- 
mined by Olofsso n et alJ Jl993h . It is also ro ughly consistent 
with the abundances of ~5 x 10~ 5 obtained bv lLindavist et al] 
(2000) for a small sample of carbon stars. Since for a given line 
intensity t he derived abundance scales inversely with the mass 
loss rate Joiofsson et aljfl993l) . this argues in favour of our 
adopted mass loss rate of 2x 10~ 7 M yr _1 , which is two orders 
of magnitude lower than the mass loss rate needed to model 
the CO emission (Lind qvist et a l ., in p rep.). Furthermore, as 
noted already by Olof sson et alJ Jl996l) . the HCN expansion 
velocity of 10.5 km s _1 is significantly lower than that of 16.5 
km s _1 measured from the CO lines. The fact that the 'at- 
tached' (present mass loss) envelope expands more slowly than 
the detached shell is even more pronounced in sources with 
old detached shells, e.g. TT Cyg Joiofsson et aIlll998Ll200ol) . 
where expansion velocities for the attached envelope as low 
as 4 km s~' has been measured. R Sc l is more reminiscent of 
U Cam JLindqvist et alJ IT99(tL 1 19991) . which has a relatively 
young detached shell, and whose attached envelope expands 
more slowly than the detached shell but not to the extent seen 
in sources with older detached shells. In the case of U Cam, 
however, HCN emission is detected from both the attached en- 
velope and detached shell, whereas HCN emission in R Scl is 
confined to the former. 

For the detached shell to be undetected in HCN, its 
HCN abundance there must be dramatically lower than in 
the attached envelope. Adopting the best fit CO model from 
Lindqvist et al. (in prep.), where the shell is located at r * 
11" and has a total mass of ~10~ 2 M , an upper limit of 
/hcn ^ 2 x 10~ 7 is obtained. The SEST / = 1 -> data 
provide the best constraint on this value since a higher HCN 
abundance in the detached shell would produce more emission 
in the line wings (|V - V+| > 11 km s -1 ) than what is ob- 
served. From the ATCA data alone, we obtain an upper limit 
of /hcn $ 1 x 10~ 6 in the detached shell. The most plausible 
explanation is the rapid photodissociation of HCN to form CN 
as the envelope expands away from the star. The difficulty in 
detecting HCN from sources with larger detached shells (e.g., 
S Se t, U Ant, and TT C yg) is consistent with this hypothe- 
sis dOlofsson et alJl996T) . A more extended distribution of CN 
as compared to HCN has been direc tly observed in circum - 
stellar en velopes around IRC+ 102 16 JDaval & Biegingl9 95). 
U Ca m ([Lindqvist et alJll996l) . and LP And (Lindqvist et alJ 
200(1 iLindqvist et alJ (12000) show how this phenomenon can 
be accounted for by a simple photodissociation model. 

6. Conclusions 

We have presented interferometric observations of HCN / = 
1 — > emission around R Scl that reveal a compact central 
source with FWHM of 1", although the data do not sample 
small enough scales to distinguish a Gaussian from, e.g., a 
disk or an optically thin sphere. We see no clear indication that 
the source structure varies with velocity, as would be expected 
from an expanding spherical shell, but such variations would 



be difficult to discern because the source is only moderately re- 
solved, and because of the hyperfine structure of the HCN line. 
Our interferometric data recover, to within the calibration un- 
certainties, all of the single-dish flux observed with the SEST. 

The compact HCN morphology and the narrow line profile 
of the higher-transition HCN lines indicate that the HCN emis- 
sion is associated with the present mass loss wind. We estimate 
a mass loss rate of 2 x 10~ 7 M yr _1 , significantly less than that 
needed to account for the detached shell inferred from single- 
dish CO data. This supports the idea that the star has undergone 
a recent (within ~10 3 yr) episode of intense mass loss. The low 
HCN abundance in the detached shell, which we estimate to be 
/hcn $2x 10~ 7 , is consistent with the rapid photodissociation 
of HCN into CN as it expands away from the star. 

Modelling the present mass loss envelope using a detailed 
radiative transfer code, the ATCA data suggest an HCN abun- 
dance and envelope size of 4.5 x 10~ 5 and 5.5 x 10 15 cm, re- 
spectively. However, such a model cannot be reconciled with 
single-dish observations of higher J transitions. The single- 
dish constraints are much better met using an HCN abundance 
of 1 .2 x 10~ 5 and an HCN envelope size of 1 .4 x 10 16 cm. The 
discrepancy may be due to some limitation in the observational 
data, or to deviations from the adopted circumstellar model 
with regards to geometry, velocity law, dumpiness, or radia- 
tive excitation. However, it is not possible within the adopted 
circumstellar model to reconcile the interferometer and single- 
dish data by just changing the mass loss rate. 

These observations demonstrate the imaging potential of 
the Australia Telescope at millimetre wavelengths, even at this 
early stage in the system development. Future high-resolution 
observations of CN and CO emission from this object, which 
are likely to reveal a detached shell morphology, are eagerly an- 
ticipated. Such observations could be conducted with the fully 
upgraded ATCA and, at higher transitions, with the soon-to-be 
completed Sub-Millimeter Array (SMA). 
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